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Abstract 

In this paper, we are interested in the time discrete approximation of E f(Xr) when X is the solution of a 
stochastic differential equation with a diffusion coefficient function of the form \x\ a . We propose a symmetrized 
version of the Euler scheme, applied to X. The symmetrized version is very easy to simulate on a computer. For 
smooth functions /, we prove the Feynman Kac representation formula u(t,x) = E t,xf{Xr), for u solving the 
associated Kolmogorov PDE and we obtain the upper-bounds on the spatial derivatives of u up to the order four. 
Then we show that the weak error of our symmetrized scheme is of order one, as for the classical Euler scheme. 
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1 Introduction 

We consider (Xt,t > 0), the R-valued process solution to the following one-dimensional Ito stochastic differential 
equation 


b(X s )ds + a f \X s \ a dW s , (1) 

Jo 

where xq and a are given constants, a > 0 and (W t , t > 0) is a one-dimensional Brownian motion defined on a 
given probability space (f2, F. P). We denote by (Ft, t > 0) the Brownian filtration. To ensure the existence of such 
process, we state the following hypotheses: 

(HO) a e [1/2,1). 

(HI) The drift function h is such that 6(0) > 0 and satisfies the Lipschitz condition 

\b(x) - b(y )| < K\x-y\, V (x,y) e R 2 . 



Under hypotheses (HO) and (HI), strong existence and uniqueness holds for equation (1). Moreover, when xq > 0 and 
6(0) > 0, the process (X t ,t > 0) is valued in [0, +oo) (see e.g. [14]). Then (X) is the unique strong solution to 


+ a [ xyw s . ( 2 ) 

Jo 

* A previous version of this paper circulated with the title: An efficient discretisation scheme for one dimensional SDEs with a diffusion coefficient 
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X t = x 0 + [ b{X s )ds 
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Simulation schemes for Equation (1) are motivated by some applications in Finance: in [ 6 ], Cox, Ingersoll and 
Ross (CIR) proposed to model the dynamics of the short term interest rate as the solution of (1) with a = 1/2 and 
b(x) = a — bx. Still to model the short term interest rate, Hull and White [13] proposed the following mean-reverting 
diffusion process 


dr t = (a(f) — b{t)r t )dt + a{t)r^dW t 

with 0 < a < 1. More recently, the stochastic-a/3p model or SAB /7-motlel have been proposed as a stochastic 
correlated volatility model for the asset price (see [ 12 ]): 

dX t = <r t Xf dWl 
d<r t = aa t dB t 


where B t = pW} + -^/(l — p 2 )W 2 , p £ [—1,1] and (W 1 , W 2 ) is a 2d-Brownian motion. 

CIR-like models arise also in fluid mechanics: in the stochastic Lagrangian modeling of turbulent flow, character¬ 
istic quantities like the instantaneous turbulent frequency (w*) are modeled by (see [9]) 

du t = —C 3 (oj t ) ( u> t - (u t )) dt - S((uit))ui t dt 

As/ Ci{u)t) 2 wtdWt 

where the ensemble average (cot) denotes here the conditional expectation with respect to the position of the underlying 
portion of fluid and S(ui) is a given function. 

In the examples above, the solution processes are all positive. In the practice, this could be an important feature 
of the model that simulation procedures have to preserve. By using the classical Euler scheme, one cannot define 
a positive approximation process. Similar situations occur when one consider discretisation scheme of a reflected 
stochastic differential equation. To maintain the approximation process in a given domain, an efficient strategy consists 
in symmetrizing the value obtained by the Euler scheme with respect to the boundary of the domain (see e.g. [4]). 
Here, our preoccupation is quite similar. We want to maintain the positive value of the approximation. In addition, we 
have to deal with a just locally Lipschitz diffusion coefficient. 

In [7], Deelstra and Delbaen prove the strong convergence of the Euler scheme apply to dX t = k (7 — X t )dt + 
g{Xt)dWt where g : R — > R + vanishes at zero and satisfies the Holder condition | g{x) — g(y)\ < b\J\x — y |. The 
Euler scheme is applied to the modified equation dXt = k (7 — Xt)dt + g{X t A sxtXs})dWt . This corresponds to a 
projection scheme. For reflected SDEs, this procedure convergences weakly with a rate ^ (see [5]). Moreover, the 
positivity of the simulated process is not guaranteed. In the particular case of the CIR processes, Alfonsi [1] proposes 
some implicit schemes, which admit analytical solutions, and derives from them a family of explicit schemes. He 
analyses their rate of convergence (in both strong and weak sense) and proves a weak rate of convergence of order 1 
and an error expansion in the power of the time-step for the explicit family. Moreover, Alfonsi provides an interesting 
numerical comparison between the Deelstra and Delbaen scheme, his schemes and the present one discussed in this 
paper, in the special case of CIR processes. 

In section 2, we construct our time discretisation scheme for (Xt,t £ [0, T]), based on the symmetrized Euler 
scheme and which can be simulated easily. We prove a theoretical rate of convergence of order one for the weak 
approximation error. We analyze separately the cases a = 1/2 and 1/2 < a < 1. The convergence results are given 
in the next section in Theorems 2.3 and 2.5 respectively. The sections 3 and 4 are devoted to the proofs in this two 
respective situations. We denote ( Xt,t £ [0,T]) the approximation process. To study the weak error E f(Xr) — 
E J(Xt), we will use the Feynman-Kac representation E f(X^_ t ) = u(t,x) where u(t,x) solves the associated 
Kolmogorov PDE. The two main ingredients of the rate of convergence analysis consist in, first obtaining the upper- 
bounds on the spatial derivatives of u(t, x) up to the order four. To our knowledge, for this kind of Cauchy problem, 
there is no generic result. The second point consists in studying the behavior of the approximation process at the 
origin. 

Let us emphasis the difference between the situations a = 1/2 and 1/2 < a < 1. The case a = 1/2 could seem 
intuitively easier as the associated infinitesimal generator has unbounded but smooth coefficients. In fact, studying the 
spatial derivative of u(t, x) with probabilistic tools, we need to impose the condition 6 ( 0 ) > tr 2 , in order to define the 
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derivative of X'f with respect to x. In addition, the analysis of the approximation process (X) at the origin shows that 
the expectation of its local time is in At b ^ a ^ a . 

In the case 1/2 < a < 1, the derivatives of the diffusion coefficient of the associated infinitesimal generator are 
degenerated functions at point zero. As we cannot hope to obtain uniform upper-bounds in x for the derivatives of 
u(t,x), we prove that the approximation process goes to a neighborhood of the origin with an exponentially small 
probability and we give upper bounds for the negative moments of the approximation process (. X ). 


2 The symmetrized Euler scheme for (1) 


For Xq > 0, let (Xt, t > 0) given by (1) or (2). For a fixed time T > 0, we define a discretisation scheme (X tk , k = 
0 by 


J Xq=Xq 

>0, 

\ Xt k+1 = 

X tk + b(X tk )At + *Xt k (Wt k+1 - w tk ) 


k = 0,..., N — 1, where N denotes the number of discretisation times tk = kAt and At > 0 is a constant time step 
such that NAt = T. 

In the sequel we will use the time continuous version (Xt, 0 < t < T) of the discrete time process, which consists 
in freezing the coefficients on each interval [tk,tk+i). 


X t 


X v{ t) + (t- v(t))b(X vW ) + aX; {t) (W t - W v(t) ) 


(4) 


where rj(s) = supj-g^ m{tk'i tk < s}. The process (X t , 0 < t < T) is valued in [0, +oo). By induction on each 
subinterval [tk,tk+i), for k = 0 to N— 1, by using the Tanaka’s formula, we can easily show that (Xt) is a continuous 
semi-martingale with a continuous local time (Lf(X)) at point 0. Indeed, for any t £ [0, T], if we set 

Z t = X^ t ) + b(X v ^)(t — r](t)) + <rX vW (W t — W v ( t ))j (5) 


then Xt = \Zt \ and 


X t = xq 


sgn (Z s )b(X v{s) )ds + a sgn (Z s )X“ {s) dW s + -L°(X), 


( 6 ) 


where sgn(x) := 1 - 2 l( x < 0 ). 

The following lemma ensures the existence of the positive moments of (Xt), starting at xq at time 0, and of (Xt), 
its associated discrete time process: 

Lemma 2.1. Assume (HO) and (HI). For any Xq > 0, for any p > 1, there exists a positive constant C, depending 
on p, but also on the parameters b( 0), K, a, a and T, such that 


E 


sup X t 
te[o,T] 


2 V 


■E 


sup X t 

iG[0,T] 


■2 P 


<C(l + x 2 0 p )- 


(7) 


In the following proof, as well as in the rest of the paper, C will denote a constant that can change from line to 
line. C could depend on the parameters of the model, but it is always independent of At. 

Proof. We prove (7) for (Xt, 0 < t < T) only, the case of (X t , 0 < t < T) could be deduced by similar arguments. 
By the Ito’s formula, and noting that for any t £ [0, T] f* X^ P 1 dL° (X) = 0, we have 

X 2 / = x 2p 


2 p 


X 2 ? 'sgn (Z s )b(X^ s) )ds 


+2po fl X] P \gn(Z s )X a v{s) dW s + a 2 p(2p - 1) f X 

Jo 


— 2 P-‘ 2 ^TF 2ol 1 

Xn( s )d s - 


( 8 ) 
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To prove (7), let’s start by showing that 


sup < C(l + x 2 0 p ). 

te[o,T] ' 


(9) 


(7) will follow from (9), (8) and the Burkholder-Davis-Gundy Inequality. Let r n be the stopping time defined by 
T n = inf{0 < s < T; X s > n}, with inf {0} = 0. Then, 


f2p 


< xl p + 2pE 


rtATn 


7 2p-l 


x/ b(x vis) )ds + a 2 p(2p-l)E 


rtATn 


x] P 2 X 2 v a {s) ds 1 . 


By using (HO), (HI) and the Young Inequality, we get 

EX 2 t P ATn <x 2p + Tb(0) 2p + (2 p - 1)E Qf 


X 2p ds 


+ 2pKE 


—2p— 1- 


X s X v ( s )ds )+<7 p(2p~- 1)E 


T 2p ~ 2 Ar 2a l 


Replacing X s by (4) in the integrals above, by using another time (HI) and the Young Inequality, we easily obtain 
that for any t € [0, T], 


?-v 2 p 


EX^ mTn <xl P +C 


(■*/.*'■< 


ir 2p 

p(s) Ar r 


ds 


where C > 0 depends on p, />(0), K, a, a and T. A discrete version of the Gronwall Lemma allows us to conclude 
that 


sup E(x 2p kAT ) <C(l+x 2 0 p ), 

rfl M V A 


k=0,...,N 


/^F 2 P\ 


for another constant C, which does not depend on n. Taking the limit n —> +oo, we get that sup fe=0 N E{X t ) < 
C( 1 + Xq P ), from which we easily deduce (9) using (4). □ 


2.1 Main results 

In addition of hypotheses (HO) and (HI), we will analyze the convergence rate of (3) under the following hypothesis: 
(H2) The drift function b{x) is a C 4 function, with bounded derivatives up to the order 4. 

2.1.1 Convergence rate when a = 1/2 

Under (HI), (X t , 0 < t < T) satisfies 

b{X s )ds + a I 

Jo 

When b(x) is of the form a — fix, with a > 0, (X t 'j is the classical CIR process used in Financial mathematics to 
model the short interest rate. When b(x) = a > 0, (Xt) is the square of a Bessel process. Here we consider a generic 
drift function b(x), with the following restriction : 

(H3) 6(0) > cr 2 . 

Remark 2.2. When xq > 0 and 6(0) > a 2 /2, by using the Feller’s test, one can show that P(to = oo) = 1 
where tq = inf{t > 0; X t = 0}. We need the stronger Hypothesis (H3) to prove that the derivative (in the sense 
of the quadratic mean) of X/ with respect to x is well defined (see Proposition 3.4 and its proof in Appendix B). In 
particular, we need to use the Lemma 3.1 which controls the inverse moments and the exponential inverse moment of 
the CIR-like process ( X t ),for some values of the parameter v = 2b ^ — 1 > 1. 


sfX~ s dW s , 0 < t < T. 


( 10 ) 
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Section 3 is devoted to the proof of the following 


Theorem 2.3. Let f be a R-valued C 4 bounded function, with bounded spatial derivatives up to the order 4. Let 
a = t; and Xq > 0. Assume (HI), (H2) and (H3). Choose At sufficiently small in (3), i.e. At <1/(2 1\) A Xq. Then 
there exists a positive constant C depending on f, b, T and Xq such that 

f b(0) \ 

At + it?; ' 


Under (H3), the global theoretical rate of convergence is of order one. When 6(0) < cr 2 , numerical tests for the 
CIR process show that the rate of convergence becomes under-linear (see [8] and the comparison of numerical schemes 
for the CIR process performed by Alfonsi in [1]). 


2.1.2 Convergence rate when 1/2 < a < 1 

Under (HI), (X t , 0 < t < T) satisfies 


X t =x o + / b(X s )ds + cr / XfdW s , 0 < t < T. 


( 11 ) 


We restrict ourselves to the case 
(H3’) x 0 > ^ At. 

Remark 2.4. When xq > 0, the Feller’s test on process ( X t ) shows that it is enough to suppose 6(0) > 0, as in (HI), 
to ensure that P(ro = oo) = 1 ,/or To = inf{t > 0; X t = 0}. 

In Section 4, we prove the following 

Theorem 2.5. Let f be a 1R- valued bounded C 4 function, with bounded spatial derivatives up to the order 4. Let 
1 < a < 1. Assume (HI), (H2) and (H3’). Choose At sufficiently small in (3), i.e. At < 1/(4A'). Then there exists a 
positive constant C depending on f, a, a, 6 , T and Xq such that 


\Ef(X T )-Ef(X T )\<C 1 + 


„9(a) 


At, 


where q(a) is a positive constant depending only on a. 


3 The case of processes with a = 1/2 

3.1 Preliminary results 

In this section (X t ) denotes the solution of (10) starting at the deterministic point xq at time 0. When we need to 
vary the deterministic initial position, we mention it explicitly by using the notation (Xf) corresponding to the unique 
strong solution of the equation 


X? = x + [ b(Xf)ds + o f y/XfdWs 
Jo Jo 


( 12 ) 
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3.1.1 On the exact process 

We have the following 


Lemma 3.1. Let us assume (HI) and (H3). We set v = —^' — 1 > 1. For any p such that 1 < p < u, for any 
t £ [0, T] and any x > 0, 

E {Xf)~ x < C( TjaT 1 and E(Xf)~ p < C{T)t~ p or E(Xf)~ p < C(T,p)x~ p . 

2 2 

Moreover for all p < v ° , 

Eexp ^i J (Xf^ds'j < C(T ) (l + aT"/ 2 ) , (13) 

where the positive constant C(T ) is a non-decreasing function ofT and does not depend on x. 

Proof. As b(x) > b(0) — Kx, The Comparison Theorem gives that, a.s. Xf > Y t x , for all t > 0, where (Y t x , t < T) 
is the CIR process solving 

Yf = x + f (6(0) — KYf) ds + a f ^/YfdW s . (14) 

Jo Jo 

In particular, Ecxp(/r f ( X x )~ 1 ds) < Eexp (p f ( Yf)~ 1 ds ). As 6(0) > er 2 by (H3), one can derive (13) from the 
Jo Jo 

Lemma A.2. Similarly, for the upper bounds on the inverse moments of Xf, we apply the Lemma A. 1. □ 


3.1.2 On the associated Kolmogorov PDE 

Proposition 3.2. Let a = 1/2. Let f be a R -valued C 4 bounded function, with bounded spatial derivatives up to the 
order 4. We consider the WL-valuedfunction defined on [0, T] x [0, +oo) by u(t, x) = E f(Xif_ t ). Under (HI), (H2) 
and (H3), u is in C 1 , 4 ([0, T] x [0, +oo)). That is, u has a first derivative in the time variable and derivatives up to 
order 4 in the space variable. Moreover, there exists a positive constant C depending on f, 6 and T such that, for all 
x £ [ 0 , +oo). 




sup 

te[o,T] 

U LL . 

<C(1 + x ), 

(15) 

4 

l u IIl°°([0,T1x[0,+oo)) + ^ 

d k u 

dx k 



<C 

(16) 

k= 1 

L°°([0,T] x [0,+oo)) 



and u(t, x) satisfies 


du, . , . . du, , cr 2 

-{t,x) + Hx)-{t,x) + - 

u(T,x) = /( x), x £ [0,+oo). 


d 2 i l 

X Q^^ X )=^ 


(t,x) £ [0, T) x [0,+oo), 


(17) 


In all what follows, we will denote || ||oo the norm on L°° spaces. Before to prove the Proposition 3.2, we introduce 
some notations and give preliminary results: for any A > 0 and any x > 0, we denote by (X[(X), 0 < t < T), the 
[ 0 , +oo)-valued process solving 


X\ 


: (A) =x + Xa 2 t+ [ b{Xf{X))ds + a [ y/X x {X)dW s 
Jo Jo 


( 18 ) 


Equation (18) has a non-exploding unique strong solution. Moreover, for all t > 0, Xf(X) > Xf. The coefficients are 
locally Lipschitz on (0, +oo), with locally Lipschitz first order derivatives. Then Xf(X) is continuously differentiable 
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(see e.g. Theorem V.39 in [16]), and if we denote Jf( A) = the process (Jf(X),0 < t < T) satisfies the 

linear equation 


J?( A) = 1 + T J s *(A)&'(JC(A))ds + f Jf( A) 

Jo Jo 


adW s 


o SV 2 V / XF(A)' 


(19) 


By Lemma 3.1, the process (/„ —y=—, 0 < i < T) is a locally square integrable martingale. Then, for all A > 0, 
Jf (A) is given by (see e.g. Theorem V.51 in [16]), 


■TT(A) = exp b'(Xg(X))ds + ^ 


dW* 


T 2 /-t 


ds 


/o \/(A) 8 7o X*(A) y 

Lemma 3.3. Assume (H3). The process A), 0 < t < T) defined by 


( 20 ) 


Mf(A) = exp (| £ 


dW. 






8 Jo ; 

is a P -martingale. Moreover, sup te [ 0 T j E ( Jf (A)) < C where the positive constant C does not depend on x . 

Proof. By Lemma 3.1, (Mf( A), 0 < t < T) satisfies the Novikov criterion. Under (H2), 1/ is a bounded function and 


E[J?(A)]=E 


ex? (yf*b'(Xs(ty)ds S J Mf{A) 


< expdl&'HooT). 


Let (Z t (A ’ A+ 2 \o <t <T) defined by 


r (A,A+^) 


= ex P -f f c 


2 Jo v^W 


dX*( A) 




-ds 


zL r* ds 

8 JO JVf(A) 


□ 


( 21 ) 


By the Girsanov Theorem, under the probability 


llA+i 


2 such that 




dp 


7 (S,S+^) 


, the process (13 




rt d _ b(X s (a))+(a+ 2 ) ct ^ [0, Tl) is a Brownian motion on (fl, J-t,© A+ ^). Indeed, we have that 

J 0 a^xf(X) cr^xfiX) ^ V ’ 

X?{\) = x + (X + \)o 2 t + [ b(Xf(X))ds + <7 [ ^/Xf(X)dBs + 1 \ 

* Jo Jo 


Hence, 


A+i 


(X X (A)) = ^(^(A+l)) and from the Lemma 3.3, zj. X ’ X+2 ' > = exp(—| J 


a rt dBg 


_ <j^_ ft ds 

2 J o y/x7(X) 8 Jo X sW 


■) 


is a 1 


tA+j 


-martingale. The following proposition allows us to compute the derivatives of u(t, x). 


Proposition 3.4. Assume (HI), (H2) and (H3). Let g{x), h{x) and k(x) be some bounded C 1 functions, with bounded 
first derivatives. For any A > 0, let v(t, x) be the WL-valued function defined, on [0, T] x R[j_, by 


}{t,x) = E 


g(Xf( A))exp(/ k(Xf(X))ds) 


E 


h{Xf(X)) exp( / k{X^X))d6) 
Jo 


ds. 


Then v(t, x) is of class C 1 with respect to x and 


f^-E 


exp (jf k(x;m)d,) L'(Xf(A))Jf (A) + j(Xf (A)) J k'(X;(\))j;(X)d, 


E 


ex p( / k(X$(\))d9) h'(Xf(X))Jf(X) + h(X*(X)) / k'(Xg(X))J%(X)d6 


ds. 


1 








































The proof is postponed in the Appendix B. 

Proof of Proposition 3.2. First, we note that u(t, x) = E f(Xif_ t ) is a continuous function in x and bounded by ||/||oo- 
Let us show that u is in C 1,4 ([0, T] x [0, +oo)). / being in C” 4 (M), by the Ito’s formula, 

n(t, x) = f{x) + £ ‘ E (b(X*)f'(X*)) ds + ^-jy E (Xff"(Xf)) ds 

+ aE^ T ‘ y/X*f’(XZ)dW, 

f is bounded and (Xf) have moments of any order. Then we obtain that 

y,x) = -E • 

Hence, ^ is a continuous function on [0, T] x [0, +oo) and (15) follows by Lemma 2.1. 

By Proposition 3.4, for x > 0, u{t, x ) = E f(X^_ t ) is differentiable and 

fioi 

— (f,x)=E(/W_ t (0))J^_ t (0)). 

Hence, by using the Lemma 3.3, x)| < H/'HooE (J^_ t (0)) < CH/'Hoo. We introduce the probability Qs suc h 

= | . Denoting by E 5 the expectation under the probability Q 5 , we have 

y, X ) = E^ (nxZ_ t (0))Z$$J$_ t (0j) . 

From (20), as Wt = B? + I * — " ds, we notice that 
* J 0 2y/xf(p) 




■Pity = exp ( / b '(X*(0)ds+- 


dB} 


ds 


and that Jf_ t {0) = exp (/ Q T t b'{Xf (0))ds), 


E 2 


r T-t 


f( X T- t ( 0))exp 


b'(X*(0))ds 


following expression for 4^(f, x): 


o 

du t 


2 Jo y/xf(fi) 8 Jo X S(0) y 

from the definition of z}°' 2 ^ in (21). Hence, ^{t,x) = 
. As PM 2 (X x (0)) = Cp{X X (B)), for x > 0, we finally obtain the 


du, . 

^((,*)=E 


p(y T V(x,‘(i))i 


( 22 ) 


Now the right-hand side of (22) is a continuous function on [0, T] x [0,+oo) so that u € C' 1 , 1 ([0,T] x [0,+oo)). 

/ 


Moreover for x > 0, by Proposition 3.4, x) is continuously differentiable and 


&£(*,*) = E 


+E 


f'\ X T-t(\W-t{\) J q V(XZ(\))ds 


(23) 
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As previously, we can conclude that 


{t, x ) 


is bounded uniformly in x. In order to obtain an expression for 


^4£(f, x) continuous in [0, T] x [0, +oo) and also to compute the third derivative, we need to transform the expression 
in (23) in order to avoid again the appearance of the derivative of ,/f (|) that we do not control. Thanks to the Markov 
property and the time homogeneity of the process (|)), for any s £ [0, T — t], 


E 


2 > 

T-t 


f'( X T-t(n)) ex P 


&W(o)) / ^ 


= E 


pT—t—s 


f'{Xj,_ t _ s (-)) exp 


b'm-)) 


V=X?(})) 


By using (22), we get E[/'(X£_ t (i)) exp( J b'(X*{^)))/F s ] = + s,X s x (i))). We introduce this last 

equality in the second term of the right-hand side of (23): 


E 


r T-t 


/'(^T-t(o)) ex P 


b'(x:u)du 



T-t 


b n mh)jf{hda 


= E 


£ * E ^/'(^_ t (i))ex P ^ T V(^(i))d«j j £ 
x exp f I* b\XZ(h)du) b"{X*{ l -))J*{ l -)ds 


pT-t 


E 


g(f + s,X*(i)))exp( / & , (X;Q)d U )6"(X“(^))^(^) 


ds. 


Coming back to (23), this leads to the following expression for a:): 


r T-t 


E 


d 2 u 

a? (f - l) = E 

du, 


pT-t 


nX°_ t (-))J*_t(-)ex p 


,i. 


&'(X*(-))ds 


— (f + s,**(-))) exp ( / b'(XZ(-))du ) b"{X*{-))J*{-) 


1 , 


ds. 


We introduce the probability Q 1 such that 


= -tE)-T hen, 


Ft z < 




(xf _,(!)) exp 


/•T-t 


b'(X*(-))ds 


pT-t 


E 1 


Ahd) du , 


— (f + s,X*(-)))exp( / b’{X*{-))du ) 6"(X*(-))J*(-) 




ds. 


Again for all d € [0,T], we have that Zg 2 ' 1 ' 1 Jg (|) = exp ^ iff b'(X*(^))duj and 


S3 (*.-)= ei 


pT—t 


f"(X£_ t (-))ex p 2/ b'(Xf( — ))ds 


+ 


rT-t 

/ E 1 

Jo 


^(f + S ,X-(i)))exp(2 / 6'(^Q)dtt)&"(X*Q) 


ds. 
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As (X x (i)) = we finally obtain the following expression for =4f (t, at): 


a 2 « 

dx 2 


(f, at) = E 


/"(X£_ t (l))exp V(X?(l))d. 


+ fo~ 


g(i + S ,X*(l)))exp (2 j\\X^{l))d^j 6"(X*( 1)) 


(24) 


ds 


from which, we deduce that u £ C' 1,2 ([0, T] x [0, +00)). As JJ(1) = exists and is given by (20), for * > 0, 


(f, x) is continuously differentiable (see again Proposition 3.4) and 


d 3 u 
dx 3 


(25) 


(t,x) = E^exp^2 J b'(X£(l))ds^ 

x [/( 3 )(Xf_,(l))Jf_ t (l) + 2/"(Xf_ t (l)) / 0 T -‘ 6"(Xf (1)) J*(l)ds] } 

+ /o T_t E {exp (^2 jT 6'(X*(l)du) 

x [ J *(!) (&(* + s, A7(1)))6"(X-(1)) + f±(f + s,X®(1)))&( 3 )(X*(1))^ 

+2§|(f + a,X; (1)))6"(X*(1)) / 0 S 6"(X*(1)) Jf(l)du] } da. 

f ) 3 

By Lemma 3.3, the derivatives of / and b being bounded up to the order 3, we get immediately that | §^§(t, x)\ < C 
uniformly in x. 

The computation of the fourth derivative uses similar arguments. We detail it in the Appendix C. 

In view of (15) and (16), one can adapt easily the proof of the Theorem 6.1 in [10] and show that u(t, x) solves the 


Cauchy problem (17). 


□ 


3.1.3 On the approximation process 

According to (3) and (6), the discrete time process (X) associated to (X) is 

Xo = x 0 , 

Xt k+1 = 


x t k + b(X tk )At + <tJX tk (W tk+1 - W tk ) 


k = 0,..., N — 1, and the time continuous version (X t , 0 < t < T) satisfies 

X t =x 0 + J sgn(Z s )b(X v ( s) )ds + a J sgn (Z S )^J X v{s) dW s + ^£?(X), 
where sgn(x) = 1 — 2 l( a .< 0 ), and for any t £ [0, T], 


Zt — X,j( t ) + (t — rj(t))b(X ^(t)) + cry X v ^(W t — W v (t))- 


(26) 


(27) 


(28) 


In this section, we are interested in the behavior of the processes (X) and (Z) visiting the point 0. The main result is 
the following 

Proposition 3.5. Let a = ^. Assume (HI). For At sufficiently small (At <1/ (2K)), there exists a constant C > 0, 
depending on 6(0), K, a, xq and T but not in At, such that 


E (l°(X) - L° (t) (X)/ X vW ) < CAtexp (- X??(t) 


and E Vi(X) < C ( — ) 
Xq J 


l 16(J 2 At ) 

Ho) 


(29) 
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The upper bounds above, for the local time (L°(X)), are based on the following technical lemmas: 

Lemma 3.6. Assume (HI). Assume also that At is sufficiently small (At <1/(2 K) A Xq). Then for any 7 > 1, there 
exists a positive constant C, depending on all the parameters 6(0), K, a, Xq, T and also on 7 , such that 


sup 

k= 0 ,...,N 


Eexp 


X tk \ 

ya 2 At.) 


< C 



Lemma 3.7. Assume (HI). For At sufficiently small (At < 1 /(2K)), for any t £ [0,T], 


P (Z t < 0/ X„ (t )) < ^ exp 


X. 


•n(A) 


2(1 - KAt)~ 2 



As 2(1 — KAt) 2 > 1 when At < l/(2/\), the combination of Lemmas 3.6 and 3.7 leads to 

P(Z t <0) <c(^j . (30) 

We give successively the proofs of Lemmas 3.6, 3.7 and Proposition 3.5. 

Proof of Lemma 3.6. First, we show that there exits a positive sequence (pj . 0 < j < N ) such that, for any k £ 

{I,---,N}, 


Eexp (^- ^ 2 At ) - exp exp(-x 0 p k ) ■ 

We set po = By (26), as —b(x) < —6(0) + Kx, for all x £ R, we have that 

Eexp (-p 0 X tk ) < Eexp (^-p 0 ^4., + (6(0) - KX th _f)At + ayjX tk _ ± AW tk ^j 


' a2 . 2 , 


A W tk and X tk _ 1 being independent, Eexp(-p 0 ayJ X tk _ 1 AW tk ) = Eexp(— p( ) AtX tk _ 1 ). Thus 

Eexp (—p 0 X tk ) < exp (-/r o 6(0)Af) Eexp (^-p 0 Xt k _ k (l-KAt- ^-poAt^j 
= exp (-p o b(0)At) Eexp (-piXt^) , 

where we set pi = po(l — KAt — ^-paAt). Consider now the sequence (/Xj)j £ N defined by 

( To = 7<T *At ’ 

[Pj = pj-i (l - KAt - ^pj-xAtj , j> 1. 


(31) 


An easy computation shows that if 7 > 1 and At < then (pf) is a positive and decreasing sequence. For any 
j £ {0 ,... ,k — 1}, by the same computation we have 

Eexp (- pjX tk _ d ) < exp (-b(0)pj At) Eexp (-p j+ iX tk _ j _ 1 ) 

and it follows by induction that 

k -1 


Eexp (-p 0 X tk ) < exp [ -6(0) ^ pjAt J exp ( -x 0 p k ). 

3 =0 


(32) 
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Now, we study the sequence (pj, 0 < j < N). For any a > 0, we consider the non-decreasing function f a (x) := 

1 ^ ax , x £ R. We note that (f a o fp)(x) = f a +p{x). Then, for any j > 1, the sequence (/r ? ) being decreasing, 
2 

Hj < Hj- 1 — and 


H < fsq. At^J-l) ^ fsq .At (f^AtiN- 2 )) <■■■< 

The next step consists in proving, by induction, the following lower bound for the Hj : 

1 \ „ / Af(j -1 )/j, 0 


Hj > hi 


- K 


, Vj > 1. 


1 + ^-At(i - i)/u 0 y \i + - i) / u 0/ 

(34) is clearly true for j = 1. Suppose (34) holds for j. By (31) and (33) 

Mj+i = - Y AtfJlj ) ~ KAt ^j - yAf/^^^o)^ - KAtf^_ jAt (no) 

1 + ^-At(j “ 1 )/ i o \ _ p ( At/ip \ 


>Fi 


1 + At j Ho 


(33) 


(34) 


1 + ZyAtjHo J 


and we conclude by using (34) for Hj■ Now, we replace /r 0 by its value l At in (34) and obtain that Hj > 
At 7 J( 27 ~i+j) - % for an y j > °- Hence > 


y A , . . J_ y 2 7 -1 _ 2Ktk J_ f k 2 7 -1 _ 2fCT 

“J ' J ~ 7(J 2 27 — 1 + j cr 2 — 7a 2 Jq 27 — 1 + u a 2 

27-I /27-l + fc\ 2 FfT 

— 7cr 2 \ 27 — 1 y cr 2 

Coming back to (32), we obtain that 

<e X p(f.(0)^)e xp (x„^) 

2b(0) /n 1 x 

/ 2 7 -l ( XQ ( 2 7 -l) \ 

\ 27 — 1 + k J y Aiycr 2 (27 — 1 + k) J 

Finally, we use the inequality x a exp(— x) < a a exp(— a), for all a > 0 and x > 0. It comes that 


E exp - 


X- 


tk 


7 a 2 Af 


2KT \ ( 2 K\ 

<exp I 6 ( 0 ) I exp I (Coy I 




2b(0) 


a-*) 


exp 


27 


<c 


(£) 


2b(0) 


(i-&) 


where we set 


26 ( 0 ) 1 , 

C = ( 6 ( 0 )( 2 7 - exp 


b(0)(KT — 1 + —) + xo-ff 


□ 
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Proof of Lemma 3.7. Under (HI), b(x) > fc(0) — Kx, for x > 0. Then, by the definition of (Z) in (28), 


1 (Zt < 0 ) < P \w t - w v{t) < 


~ x v(t ){ 1 ^ K{t - nit))) - 6 ( 0 )(f - ■qit)) 


oJX 


X v (t) > o | . 


vW 


By using the Gaussian inequality P(G < (3) < 1/2 exp(—/3 2 /2), for a standard Normal r.v. G and /3 < 0, we get 


’ {Z t < 0) < -E 


exp - 


(Xn( t ){l -K{t- n{t))) + b( 0 )it - 77(f))) 2 


2a 2 it-n{t))X ri ( t ) 


1 {X, (t) >0} 


1 X 

from which we finally obtain that P (Z t < 0) < -E[exp(--- „ „ 

J v - ’ ~ 2 L 2(1 - KAt)~ 2 a 2 At 


)]■ 


□ 


Proof of Proposition 3.5. By the occupation time formula, for any f £ [tk, tk+\ ), 0 < k < N, for any bounded 
Borel-measurable function </>, P a.s 

L m mm - Lim) d Z = [ m/z) - Lum d Z 

JR JR 

= f f(Z s )d(Z 1 Z) s =a 2 f f(Z s )X tk ds. 

Jt k Jtk 

Hence, for any x > 0, an easy computation shows that 

L mm ~ Li k m/ {x tk = x }) d Z = ^ f xe mz s )/ {x tk = *}) d S 

JR Jt k 


<t>{z) 


yfx 


t k v / 27r(s — t k ) 


exp 


iz — x — bix)is - t k )) 2 


2 er 2 x(s — tk) 


ds dz. 


t k y/2itis-t k ) 

At y/E 


Then, for any tel, 

e mm - Lim/ [x tk = x })=a f 

Ju 

In particular for Z = 0 and f = tk + \, 

E { L l +1 m - L im/ =*}) = * l 

From (HI), 6 (x) > —Kx, with K > 0. Then, 

E(L° tk+i m-L 0 tk (X)/T tk )<oj o 

For At sufficiently small, 1 — KAt > 1/2 and 


y/x ^ f iz - X - 6 (x)(s - ffc )) 2 


2 a 2 x (s — ffc) 


ds. 


A* y X 


y/2ns 


y[2/KS 


exp 


exp - 


(x + &(x)s ) 2 
2 tr 2 x s 


ds. 


Xt k jl -Ksf 
2o 2 s 


ds. 


E 


(Ll +i iX)-L° tk iX)/x tk )<aJ o 


At y X 


y/2irs 


exp 




t k 


8 er 2 Af 


Now we use the upper-bound aexp(—^-) < 1, Va £ R, to get 


E 


fijx) - Lim) 


) < o - 2 

r 2 /? e [ 

/ - 

Jo V™ 


exp 


4cr 2 Af . 

< -=— sup Eexp I — 


X 


tk 


V n k=0,...,N 

We sum over k and apply the Lemma 3.6 to end the proof. 


16er 2 Af J 
Xtk 


ds. 


ds 


16cr 2 Af J 


□ 
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3.2 Proof of Theorem 2.3 


We are now in position to prove the Theorem 2.3. To study the weak error E f{Xr) — E f(Xr), we use the Feynman- 
Kac representation of the solution u(t, x) to the Cauchy problem (17), studied in the Proposition 3.2: for all (t, x) G 
[0, T] x (0, -poo), E/(Xf._ t ) = u(t, x). Thus the weak error becomes 


E f(X T ) - E/(I T ) = E (u(0, x 0 ) - u(T, X T )) 

with (X) satisfying (27). Applying the Ito’s formula a first time, we obtain that 


E [m(0,x 0 ) - u(T,X t )\ 


du , — 


dt 


du, — 


— d 2 u, — 


(s> X a ) + sgn (Z s )b(X v ^) — (s, X a ) + —A,X a ) 


-f 

-E J S gn(Z s )ayfX^^(s,X s )dW s -Ej ±g( s,X s )dL°(X) s . 


ds 


j Q'U, 

From Proposition 3.2 and Lemma 2.1, we easily check that the stochastic integral (f Q sgn (Z S )-\JX v ^—(s, X s )dW s ) 
is a martingale. Furthermore, we use the Cauchy problem (17) to get 


E [u(0, xo) - u(T , Xt)] 

r T 


,— , x du , — 


— , d 2 u , — 


(P(X vW ) - b(X,)) q^(s,X s ) + — (. X vW - X a ) g^(a,X a ) 


ds 


r i f)'ij _ _ r f _ _ \ 

E J 0 2d^ {s ’ Xs)dL{x)s + J 0 ™{i<z.< 0} K x M)^x.))ds. 


From Proposition 3.5, 


E 



du 

dx 


(s, X s )dL°(X) s 


< 

du 


dx 


E (L° t {X)) < C 



On the other hand, by Lemma 3.7 for any s G [0, T\, 


_ Qn _ 

2E ( ^{z s <o}b{X v ^) — (s, X s 


< 


du 

dx 


E 


(6(0) + KX v ^)exp 


Xy(s) 

8a 2 At 


As for any x > 0, xexp(— 16 c 3 At ) < 16cr 2 Af, we conclude, by Lemma 3.6, that 



_ _ 

^-(Z s < 0 }b(Xri(s))-g^{ s ,X s ) J ds 



Hence, 


|E [u(0, Xo) — u(T , At)] | 


< 


E 


(b(X v{s) ) - b(X s )) q^(s,X s ) + ^ (A„ w - Xg) ^(s,X s ) 


dx 2 


ds 


C 


/At 

V^o 


6 ( 0 ) 
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By applying the Ito’s formula a second time (jf is a C 3 function with bounded derivatives), 


E 

= E 


rjii 

(b(X s )-b(X v{s) )—(s,X s ) 


/ Sgn (Z e )b(X^ a )) 

Jri(s) 


fftl! _ f)l I 

(b(X e ) - b(X v{s) ) —{s , X,) + b'(Xg)—(s, x e ) 


d6 


' ® [ o X v (s) 

Jri(s) * 


_ _ fft 7 / _ cftii _ r)n 

(b(X e ) - b(X v{s) ) — (s, X 6 ) + 2 b\X e ) — {s, X e ) + b"(X e )—(s, X e ) 


dO 


+ E 


1 


h(s) 


_ , d 2 


du, 


(6(0)-6(X 1?(s) )^( s ,0)+6'(0)-( S ,0) 


dx 


dLg(X) 


so that 


E 


_ _ f)n I _ 

(b(X s )-b(X v{s) )—(s,X s ) 


cAt ( i + sup e|a' 9 | 2 + e{(i + |a: i , w |)(l2(a:)-l° w (x))} 


< 

o<s<t 

and we conclude by Lemma 2.1 and Proposition 3.5 that 


E 


_ _ f)ll _ 

( b(X s ) — b(X v ^ s ' ) ) —(s,X s ) 


< C | At + [ — 

X 0 


6 ( 0 ) 


By similar arguments, we show that 


E 


_ _ (fiii 

(X s — X v ^) X B ) 


<C \ At + 


At 

x 0 


6 ( 0 ) 


which ends the proof of Theorem 2.3. 


4 The case of processes with 1/2 < a < 1 


4.1 Preliminary results 

In this section, (X t ) denotes the solution of (11) starting at xq at time 0 and (X[), starting at x > 0 at time 0, is the 
unique strong solution to 


Xf = x + / b{X x s )ds + (7 / (X x ) a dW s 


(35) 


4.1.1 On the exact solution 

We give some upper-bounds on inverse moments and exponential moments of (X t ). 

Lemma 4.1. Assume (HI). Let x > 0. For any 1/2 < a < 1, for any p > 0, there exists a positive constant C, 
depending on the parameters of the model (35) and on p such that 


sup E 
te[o,T] 




< C{1 + x~ p ). 
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Proof. Let r n be the stopping time defined by r n = inf{0 < s < T; Xf < 1 /n}. By the Ito’s formula. 


E[(Xf AT J-P] =x~P~pE 


f* ATn b(Xf)ds 

to WF 


^ +rK J. e V(x,- a ,j» 


+ p(p + 1 ) —E 


ds 


ds 


•E 


Uo 


p(p + 1 )' 


2 (X^)p+ 2 ( 1_ “) P (Xf)P +1 


J 0 (X)f) p+2(1_a) _ 

6 ( 0 ) 


ds 


It is possible to find a positive constant C such that, for any x > 0, 


p(p + 1 )' 


‘ -,«< c 


2 a;P+2(l-a) ^P+l / 


An easy computation shows that C_ = p (2a — 1) %- (p + 2(1 — a)) ^qJ 
upper-bound above. Hence, 


p+2(l —cO 
2 q —1 


is the smallest one satisfying the 


E[(Xf AT J-P] <x-r + CT+pK f sup E [(X£ At J“ p ] ds 

Jo ee[o,s] 


and by the Gronwall Lemma 


sup E [[Xf ATn )- v ] < ( x-r + CT) exp(pKT). 
te[o,T] 

We end the proof, by taking the limit n +oo. □ 

Lemma 4.2. Assume (HI). 

(i) For any a > 0, for all 0 < t < T, a.s. (Xf) 2 *- 1 ” 0 ) > Tt(a), where (r t (a), 0 < t < T) is the solution of the CIR 
Equation: 

r t (a) = x 2< ' 1 ~ a ' 1 + f (a —\(a)r s (a))ds + 2<j(l —a) f s/r s (a)d,W s 

Jo Jo 


with 


X(a) = 2(1 — a)K + 


( 2 a — 1 ) “ (a + cr 2 (l — a)(2a — 1 )) 
6 ( 0 ) 2 ( 1 -«) 


(36) 


(ii) For all p > 0, there exists a constant C(T , p) > 0 with a non decreasing dependency on T and p, depending 
also on K, 6(0), a, a and x such that 


Eexp ("i (xjyW)s c P»- 


(37) 


(iii) The process (Mf , 0 < t < T) defined by 


Mf = exp ( ao 


dW s 


T 2 r t 


ds 




2 Jo (X-) 2 ( 1 -“) 


(38) 


16 
























is a martingale. Moreover for all p > 1, there exists a positive constant C(T,p) depending also on b( 0), cr and a, 
such that 


E Uy Mf, ') ic(T ' p) ( i+ ^) 


(39) 


^te[o, T] 

Proof. Let Z t = (Xf ) 2 ( 1_ “). By the Ito’s formula, 

Z t = x 2 ^~^ + f p(Z s )ds + 2(1 - a)a [ sJY s dW s , 

Jo Jo 

where, for all x > 0 , the drift coefficient f}(x) is defined by 

/3(x) = 2(1 — a)b(x)x~ — cr 2 (l — ct) ( 2 a — 1 ). 

From (HI), b(x) > 6(0) — Kx and for all x > 0, f3(x) > j3(x), where we set 

fd{x) = 2(1 — a)6(0)x _ 2 (i- Q ) — 2(1 — a)Kx — cr 2 (l — a) (2a — 1). 

For all a > 0 and A (a) given by (36), we consider f(x) = (3(x) — a + A (a)x. An easy computation shows that f(x) 

, . Of 1 


has one minimum at the point x * = ( a(o)- 2 (i-q)a' ) 


-l) \ 2 ( 1 “ a ) 


. Moreover, 


/(**) = { 2a-l)^ (A(a) " 2(1 " - (« + ct2 (! - «)( 2a - !)) 

and when A(a) is given by (36), f(x*) = 0. We conclude that /3(x) > a — X(a)x and (i) holds by the Comparison 
Theorem for the solutions of one-dimensional SDE. As a consequence. 


Eexp 


id 


ds 


y j 0 ( a -) 2 < i -) 


< E ( exp ( p 


ds 

r s (a) 


We want to apply the Lemma A.2, on the exponential moments of the CIR process. To this end, we must choose the 
constant a such that a > 4(1 — a) 2 a 2 and p < v Wf 1 " 1 *) 4ct .for v(a) as in Lemma A.2. An easy computation shows 
that a = 4(1 — a) 2 cr 2 V (2(1 — off 1 a 2 + (1 — a)a2y/2ji) is convenient and ( ii ) follows by applying the Lemma A.2 
to the process (r t (a), 0 < t < T). Thanks to (ii). the Novikov criteria applied to Mf is clearly satisfied. Moreover, 
by the integration by parts formula. 


Mf = 


X ? 


exp 


KXf) 

1 X x 


+ a( 1 — a)■ 


2 (X^) 2 ( 1_a ) 


ds 


< ( —— ) exp(iLT) exp 


6 ( 0 ) , . a 2 1 

~ a ~xj + Q ^ (Xf) 2 f- a ) 


ds 


To end the proof, notice that it is possible to find a positive constant A such that, for any x > 0, — a a - x2 f_ a) < 

A. An easy computation shows that 



1 ) 


(1 — a ) 3 2 a cr 2 
6 ( 0 ) 2 ( 1 "“) 



is convenient. Thus, Mf < 



exp ((K + A)T) and we conclude by using the Lemma 2.1. 


□ 
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4.1.2 On the associated Kolmogorov PDE 

Proposition 4.3. Let 1/2 < a < 1. Let f be a R-valued C 4 bounded function, with bounded spatial derivatives up 
to the order 4. We consider the WL-valuedfunction defined on [0, T\ X [0, +oo) by u(t, x ) = E f(Xif_ t ). Then under 
(HI) and (H2), u is in C' 1 , 4 ([0, T] x (0, +oo)) and there exists a positive constant C depending on f, b and T such 
that 


ll M lli=“([0,T]x[0,+oo)) 


du 

dx 


< C 

L“([0,T]x[0,+oo)) 


and for all x > 0, 


sup 

te[o,T] 


du . 


< C(l + x 2a ), 


4 


and 


sup 

te[o,T] 


E 


d k u 

dx k 


(t,x) < C 



where the constant q(a) > 0 depends only on a. Moreover, u(t, x ) satisfies 


rhi (7 ^ cP 1 ii 

+ -x>^{t,x) = 0, (t,x) G [0, T] x (0, +oo), 
u(T, x) = f(x), x G [0, +oo). 


(40) 


The following Proposition 4.4 allows us to compute the derivatives of u(t, x). Equation (35) has locally Lipschitz 
coefficients on (0, +oo), with locally Lipschitz first order derivatives. Then Xf is continuously differentiable and if 
we denote Jf = '—fr, the process (.//. 0 < t < T) satisfies the linear equation 


Jt= 1+/ J x s b\Xf)ds+ aoj; 


dW„ 


W) 1 " 


(41) 


Proposition 4.4. Assume (HI) and (H2). Let g(x), h(x) and k(x) be some C 1 functions on (0, +oo) such that, there 
exist pi > 0 and p 2 > 0, 

Vx > 0, \g{x)\ + |g'(x)| + \h(x)\ + \h'(x)\ + \k'(x)\ < C (l + x Pl + -^) , 

jfc(x)| < C (l + j.2( 1 — a ) ) ■ 

Let v be the M.-valuedfunction defined on [0, T\ X (0, +oo) by 


g(Xt) exp( / k(Xf)ds) 
Jo 


v(t, x) = E 

Then v(t, x) is of class C 1 with respect to x and 


+ / E 


h(X x ) exp( / k(X§)d0) 


ds. 


I,,,.) = E 


exp (/ k(x;)d,)(g'(x;)j; + g(xr)j k'(X*)J*is 


E 


exp( / k(X§)d0) h'{Xf)Jf + h(Xf) / k'(Xg)J$d6 


ds. 


The proof is postponed in the Appendix B. 


Proof of Proposition 4.3. Many arguments are similar to those of the proof of Proposition 3.2. Here, we restrict our 
attention on the main difficulty which consists in obtaining the upper bounds for the spatial derivatives of u(t, x) up to 


18 























the order 4. By Lemma 4.1 , for x > 0, (/ Q * , 0 < t < T) is a locally square integrable martingale. Then .J'f is 

given by 


= exp (/ b'{X x )ds + aa J 


dW s 


a 2 a 2 f 4 


ds 


(1 —q) 


2 J 0 (XJ)2(1 

equivalently Jf = exp ^ f* bJX^ds'j M t , where (M t ) is the martingale defined in (38) and satisfying (39). Thus 
, we have Jf = exp ^ f* b'(X x )ds^j M t . b' being bounded, E J x < exp(iTT) and for all p > 1, 


Or 


E sup (Jf) p < C(T) 1 + 


X“P 


yte[o,T] 

By Proposition 4.4, u(t, x) is differentiable and 

— (i,x)=E[/'(X£_ t )Jf_ t ], 
Then, | x)| < H/'Hoo exp(A'T). The integration by parts formula gives 


T x — 


(x?Y 


exp (ii 


ab(X x ) a 2 a(l — a) 


X x 


2 (X“) 2 ( 1 -“) 


ds 


We apply again the Proposition 4.4 to compute (f, x): for any x > 0, 


dx 


aJf 


a(Ji 


x\2 


XI 


+,/f QT (Vra - 


ab'(X x ) ab(X x ) a 2 a(l — a) 2 


A 2 


and 


dx 2 


a du, 


(*,x) =E[/"(Xf_ t )(JfL t ) 2 ] - - — (f,x) + aE 


W ) 2 

(•/fU) 2 


X 


(*?) 


/W- t ) 


3—2a 


Jfda 


+E 


T—t 


nx x _ t )j*_ t £ '(bf’m 


ab'(X x ) ab(X x ) ct 2 q;( 1 — a ) 5 


X a 


(X -) 2 (X ? ) 


3—2a 


Jfrfs 


(42) 


(43) 


By using the Cauchy-Schwarz Inequality with Lemma 4.1 and estimate (42), the second term on the right-hand side is 
bounded by 


ll/'llooE 

te[o,T] 

< C(T) (1 + 


sup ( j n 2 f 

£[0,T] JO 


T - b " ix;) _ 
\ S J 


ab(X x ) a 2 a(l — a) 2 

w 


(*?) 


3—2a 


ds 


1 




2.2(l+a) 

By using similar arguments, it comes that 


d 2 u 

dJ? 


{t,x) 


< C{T) 1 + 


r 2+2a 


We apply again the Proposition 4.4 to compute §4r(i, x) from (43) and next (t. x), the main difficulty being the 


number of terms to write. In view of the expression of each term can be bounded by C(T)(1 + x 
where n is the derivation order, by using the Cauchy-Schwarz Inequality and the upper-bounds E sup (dj) 1 ' < 

te[o,T] 

C(T){ 1 + x~ ap ) and su Ptg[0!T] E{X?)~ P < <7(1 + x~ p ). □ 


r —2(n— 1)—na\ 
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4.1.3 On the approximation process 


When 1/2 < a < 1, according to (3) and ( 6 ), the discrete time process (X) associated to (X) is 

k = 0,..., N — 1, 


/ Xq — Xo 


{ Xt> +1 = 

X tk + b(X tk )At + aX“ (W tfc+1 - W tk ) 


Its time continuous version (X t , 0 < t < T ) satisfies 

r t 


1 


X t = x 0 + / sgn(^ s )6(X J? ( s ))ds + cr sg^Z s )X v{s) dW s + -T/(X), 
where for any £ G [0, T], we set 

= X„ (t) + (t - t 7 (f))6(X 7)W ) + aX* {t) (W t - W v(t) ), 


(44) 

(45) 


so that, for all t £ [0, T], X t = \Zt\. 

In the sequel, we will use the following notation: 

0 BX p(At) = C(T)exp(-^ T ) , 

where the positive constants C and C(T) are independent of At but can depend on a, a and 5(0). C[T ) is non¬ 
decreasing in T. The quantity (9 eX p(At) decreases exponentially fast with At. 

In this section, we are interested in the behavior of the processes (X) or (Z) near 0. We work under the hypothesis 
(H3’): xq > We introduce the stopping time r defined by 


t = inf 


s> 0;X S < ^At 


(46) 


Under (H3’), we are able to control probabilities like P(r < T). This is an important difference with the case a = 1/2. 
Lemma 4.5. Assume (HI), (H2) and (H3’). Then 


P(r<T)<0 e xp(Af). 

Proof. The first step of the proof consists in obtaining the following estimate: 

Vfc € {0,..., X}, P ( x tk < Ma(| < Oexp(At). 


(47) 


(48) 


Indeed, as b{x) > 6(0) — Kx for x > 0, for k > 1, 

-X tfe _ 1 (l-XAf)-6(0)(l-^)Af 


x *‘ A 7T Af 


< P \ w tk - W tk _, < 


aX 


>0 ■ 


tk- 


As At is sufficiently small, by using the Gaussian inequality P(G < j3) < 1/2 exp(—/? 2 /2), for a standard Normal 
r.v. G and /3 < 0, we get 

6 ( 0 ) 


Xt k < -b=rAt 


< E 


< E 


exp 


V2_. 

( (x t> _ 1 (l-lfAt) + 6(0)(l-j-)Af)‘ 


^ 2 xZA t 


1 {A' tfc _ 1 >0} 


exp — ■ 


• 2 ( 1 — 0 ;) ' 
tk— 1 

8a 2 At 


X 


b (m-j= 2 r 


exp I- ._ 2 '7 I 1/v. 


2ct 2 X: 


■2a —1 I i {X tfc _ 1 >0} 


^fc — 1 
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By separating the events ^X tk _ 1 > y/ At| and l < v/A£ j in the expectation above, we obtain 


X tk < I < exp ( - 


V2 

Now we prove (47). Notice that 

P(r<T)< Y, 

k—0 


1 


8 ct 2 A t° 


+ exp - 


mo-js y 

2CT 2 (Ai)“-5 


= Oexp(At). 


N—l 


inf A s < , X tfc > b M At ) . 

(-1 A Z 


For each k £ {0,1,..., N — 1}, by using (48) and b{x) < 6(0) — Kx , we have 

inf Z s < b M At , > ^Af 

tk<s<tk+i Z Z 

inf < b M At , x tk > x tk < b M At 

t k <s<t k+ 1 z z y2 

inf Z s < Ma*, X tfc > ^At 

t k <s<t k+ 1 2 V2 


+ J 


< P 


inf Z s < ^Ai,X tfc > 

t fc <s<t fe+ i 2 V2 


< Oexp(Af) +E<^ 1/ 


l /_ b(0) , \ P inf - 

( At fc >_ 7T At ) \0<s<At cr 


x 1 “ 6(0) — Kx 


s B s £ 


6(0)At\ 


ax'* 2ax a J 

where ( B t ) denotes a Brownian motion independent of ( W t ). The proof is ended if we show that 


x=Xt 


. a; 1 a (6(0) — Kx) b(0)At\ 

P | inf-h ---- -s + B s < 2 ^ a j < Oexp(Af), for x > —^At. 


0<s<At (7 


CTX 


6 ( 0 ) 

V2‘ 


We use the following formula (see [2]): if (B^ , 0 < t < T) denotes a Brownian motion with drift y, starting at yo, 
then for all y < yo. 


inf B£ < y 

o <s<t s ~ y 


= A* (mzx + A!') +1 


V2t y/2 J 2 


exp(2/x(y- y 0 ))erfc 


yo ~ V _ jWt\ 
V2i y/2 ) ’ 


where 


erfc(z) = -^1 exp dy, for all z£l. We set y = 


_ (b(O)-Kx) = 

rrr Q: 9 i/U 


and we choose ?/ = 


satisfying y < y 0 , if x > ^=-A t. Then 

x 1_a (6(0) — Kx) 


P inf 


o<s<A t d ax° 

M 

2 


s +VF S < 


6(0)Af\ 


2ax° 




- ierfc ( x ~^ i , ( ft (°) “ Kx ) \ 


ax' 


l V / 2Al 


ax' 


'V2 J 


1 

+ - exp I - 


2(6(0) — Kx) f 6(0)AA 
2 ~ 


b{0)At 


erfc 


ax 


' t y/2At 


(6(0) — Kx) v/A i\ 
ax a V2 ) 


:= A(x) + B(x). 

For any z > 0, erfc(^) < exp(—z 2 ). Then, if At is sufficiently small, we have 


A(x) < exp — 


(x{l~ KAt) + ^At) 2 ' 


2a 2 x 2a At 


< exp I — 


x 2(l-a) 

8a 2 At 
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and, for x > ^=-At, A(x) < exp ) = O exp (At). Now we consider B{x). If x > , then as 

for A(x), we have 


B(x) < exp ( — 


2(6(0) - Kx) 


6(0)Af \ 


)) exp (‘ 


(x - ^2^ - (6(0) - Kx)Xt) 2 ' 


a 2 x 2a 2At 


( (x(l — KAt) + ^MA£) 2 \ 

= eX P- _9-9„nA. - S eXP 


x 2(1_a) 
8a 2 At 


and B{x) < exp (-) = O exp {At,), forx > 


if 6(0) A* / T ^ |6(0)At , 2(b(0)-Kx) , 6(0)Ai\ ^ h (0) 2 At( ^ 2 ) 

11 7T At - x < (i+A'At)’ men —fa; —J A -^23- 


and 


B(x) < exp ( — 


2(6(0) - Kx) 


T 2 r r 2cx. 


x -m*) 


< exp — 


6(0) 2 Af(-A - 1) 


For x > ^j=-A t, we get B(x) < exp ( — 


2 2 “b(0) 2(1 -“ ) (Aj-i)(l+AAt) 2 

32 ^^ 2 (At) 2 S=i 


= Oexp(At). 


□ 


Lemma 4.6. Assume (HI), (H2) one/ (H3’). Let r 6e f/ze stopping time defined in (46). For all p > 0, f/jere exists a 
positive constant C depending on 6(0), a, a, T and p but not on At, such that 


Vf <£ [0, T], E 


1 


J tAr 


< c ( i + — ). 

X n 


Proof. First, we prove that 


(49) 


Vf e [o,T], 


X 

Zt < 


v(t) 


< Ocxp(Af). 


(50) 


Indeed, while proceeding as in the proof of Lemma 4.5, we have 

p (z < Iah) < Eexp ( ^r,Wa-mt-v(t))) + 2 bmt-vm 
V 2 / \ 8a 2 {t-p(t))X 2 v a (t) 

By using (a + 6) 2 > a 2 + 2a6, with a = X r/ ^ t ' ) (l — 2K(t — 77 (f))) and 6 = 26(0)(f — 77 (f)), 


Zt < 


A 




< E exp — 


X5t, a) (l - 2KAt) 2 


L ^(t) 


8cr 2 Af 


exp -- 


6(0) (1 - 2KAt) 


2cr 2 X 


v(t) 


For Af sufficiently small. 


Z t < 


X 


v(f) 


< E exp 


x n{t) \ 
32a 2 At J 


exp -- 


6 ( 0 ) 


4ct 2 A 


vW 


By separating the events j X v ^ > \fX i j and | X v ^ < y / Af } in the expectation above, we obtain 


X 

Zt < 


v (t) 


< exp ( — 


32cr 2 (Af) c 


+ exp - 


6 ( 0 ) 


4cr 2 (Af)“-2 


r = Oexp(Af). 
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Now we prove (49). Notice that Z tAr = X t/sT , by the Ito’s formula 


j-iAr X 

—— = — -pl '—v+i' ds-pcr / 

Z P tAT 4 J o Z p s +1 Jo Z-" 

Taking the expectation and using again b(x) > 6(0) — Kx, we have 


1 1 f b(X {s) ) r” X r,(e) „„ , , ^y[™X. 

— = ^p~P +1 ds-pa dW s +p(p +1)— — 

tA r x Q j o % j 0 2. / Jo Z 


ri(s) 


■P+2 


ds. 


E 


1 


Z 


tAr 


<4r ~P E 


rtAr 


6 ( 0 ) 


rtAr 


x 


^ ds ) + ” KR [i # 


ds 


2 ( r 

+p(p+ i )^- e ( 


tAr x 2 ° 


p(s) 


z 


■P+2 


ds . 


By the definition of r in (46), 

/ /*£Ar “F 


E 


l p( 


< 


\d 0 Z 

/ 2 y +1 


, ds ) = IE / 1 x , 

:p+1 I [Jo (z s <^)z p+1 


^>ds +E f T i x () ^rds 

/ Wo (^>y^)z p+1 , 


) T sup 
/ te[o,T] 


M 1 / 2 


Zf < 




v(t) 


\6(0)Af 

We conclude, by the Lemma 2.1 and the upper-bound (50) that 

/ rtArx 


sup 

tG[0,T] L 


E 


(•*4>) 


E 


/ o Z p+1 


7?(s) ds) <C + 2 I E ( W_ 


1/2 + 2 [ e(W— ) ds. 
do \Z sAt 


ds. 


Similarly, 


E 


rtAr x 2a 


rf(s) 


—2a 


/ _ . „ ds = E / 1 _ x , - 

I o Z P / Wo (^<^)y +2 


-wr* 


—2a 
A p(s) 
(Z B >^1) 2 P+2 


ds 


( ptAr 

l 1 


x 


p( s ) 


rtAr 




ds + 2 2 “E 


ds 


-p+2(l-a) 


By using again the Lemma 2.1 and the upper-bound (50), we have 

P+2 


E 


Finally, 


1 


—2a 

A p(+ 


(z„<^l)2 P+2 


ds <T 


< T 


6(0) Ay te S [0 P T] 

2 \ p+2 


\ 


2 ) 

(&(o)a t) 0exp(At) - a 


Z t <^#l 1,/E 


(*?«>) 


E - 


rtAr 


- r? + E 


p6(0) 2 2 “ 1 p(p+ 1)(T 2 


zW 1 


—p+2(l-a) 


ds + C E 




ds + C. 


We can easily check that there exists a positive constant C such that, for all z > 0, p ° 

Hence 


E 


Z 


l \ l r 

<—+C I E 


J tAr / 0 

and we conclude by applying the Gronwall Lemma. 


1 


Zj „ 


ds + C 


< C. 


□ 
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4.2 Proof of Theorem 2.5 


As in the proof of Theorem 2.3, we use the Feynman-Kac representation of the solution of the Cauchy problem (40), 
studied in the Proposition 4.3: for all (t, x) £ [0, T\ x (0, +oo), E f(X^_ t ) = u(t, x). Thus, the weak error becomes 


E f(X T ) - E f(X T ) = E (u(0, x 0 ) - u(T, X T )). 


Let r be the stopping time defined in (46). By Lemma 4.5, 


E(u(TAt,Z TAt )-u(T,X t )) 

— ~ IMIz,<*>([0,T]x[0,+oo]) ^ (r < T) < ^exp (At). 

We bound the error by 

|E (u(T, X T ) - u(0, x 0 )) | < |E (u(T A r, Ztat) ~ u( 0, x 0 )) \ + O exp (At) 

and we are now interested in E ( u(T A r, Ztat) — (it(0, Xo)). Let L and C z the second order differential operators 
defined for any C 2 function <j(x) by 


= 6(x)^(x) + yx 2 “|-|(x) and C z g(x) = b(z)^-{x) + ^-z 2a ^(x). 

From Proposition 4.3, u is in C 1,4 ([0, T] x (0, +oo)) satisfying x) + Lu(t, x) = 0. X t has bounded moments 
and the stopped process (XtAT)=(ZtAr) has negative moments. Hence, applying the Ito’s formula, 

pT At 

E' — 


[u(T A r, Xtat) - m(0,x 0 )] = E J [c^^u - Luj (s, X S )ds, 


Notice that 


MC z u - Lu) + b( z )d x (C z u Lu) + °-z^dUC z u - Lu) 
= C 2 z u — 2C z Lu + L 2 u 


and by applying again the Ito’s formula between r)(s ) and s to ^ Cjf ( } u — Lu^j (s, X s ), 


E [u(T A r, X T at) ~ w(0, x 0 )] 

i-T PS 


E 


10 Jr)(s) 


l(e< T ) f Cy t u — _Lu + L 2 u 


L n(s) 


Tj(s) 


) (9,Xg) 


d6ds. 


(. C 2 u — 2 C z Lu + L 2 u) ( 9, x) combines the derivatives of u up to the order four with b and its derivatives up to the 
order two and some power functions like the z 4a or x 2a-2 . When we value this expression at the point (z. x) = 
(X v (sAt)i Xg A r), with the upper bounds on the derivatives of u given in the Proposition 4.3 and the positive and 
negative moments of X given in the Lemmas 2.1 and 4.6, we get 


E 


h 


L ^ 

which implies the result of Theorem 2.5. 


u — Lu + L z u 

Xri(s) 


) (0,Xg) 


< C 1 + 


M a ) 


A On the Cox-Ingersoll-Ross model 

In [6], Cox, Ingersoll and Ross proposed to model the dynamics of the short term interest rate as the solution of the 
following stochastic differential equation 

f dr* = (a — brf)dt + a^/rfdWt, 

\ rl =x > 0, y ’ 

where (W t . 0 < t < T) is a one-dimensional Brownian motion on a probability space (12, F. P), a and a are positive 
constants and b £ R. For any t £ [0, T], let Ft = cr(s < t, W s ). 
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Lemma A.l. For any x > 0 and any p > 0, 


E 


( r x)p 


i ( 2b y 

T(p) y <r 2 (1 —e~ 6t ) ) 


X fg (2 - 0)^-!(l - exp dd , 


where T(p) = / 0 + °° 1 exp(—p > 0, denotes the Gamma function. Moreover, if a > <r 2 


E 


Q bt 


< 


and, for any p such that 1 < p < =§ — 1 

1 


E 


t\ p 


1 ( 2e |b|: , 

< _. . I —-— I orE 


_(rf)pj T(p) V ° 2 t 




<C(P,T)-, 


(52) 


(53) 


(54) 


where C(p , T) is a positive constant depending on p and T. 

Proof By the definition of the Gamma function, for all x > 0 and p > 0, x~ p = L(p) -1 / 0 + °° m p_1 exp (-ux)du, so 
that 


E 


1 


F+OO 


_(rf)Pj T(p) J o 
The Laplace transform of rf is given by 

Eexp(— urf) = 


u p x Eexp (-ur^)du. 


uL(t)((t, x) 


(2uL(t) + l) 2a /°- 


■ exp - 


2 uL{t) + 1 ) ’ 


where L(f) = ^(1 - exp(-6f)) and ({t,x) = a 2 {e J^ lt )_ 1) = xe bt /L(t), (see e.g. [15]). Hence, 


E 


1 

1 

[+°° uP- 1 

F wL(t)C(f,x)\ 

l(r?) p \ 

r (p)J 

0 (2 uL{t) + l)2a/a 2 C “ P 

i, 2uL(t) + 1 y 


du. 


By changing the variable 6 = 2 2 uL{t)+i 


in the integral above, we obtain 


E 


1 


l(r?) p \ 2 p T{p)L(t) 
from which we deduce (52). Now if a > o 2 , we have forp = 1 


— f 1 (2 — 9)9 P ~ 1 (1 — 8)^~ p ~ 1 exp ( 

0L(f)P J 0 y ’ ^ ’ P \ 2L(t) J 


dd, 


and for 1 < p < |f — 1, E 

(1_ e -|«>|t) > Ifolte-I 6 !*. 


1 


< 


1 


< 


2L(t) 
1 


exp 


n —bt 


2L{t) 
2 p \b\ p 


d0< — 
x 


_(rf) p J “ 2 Pr(p)L(f)P cr 2 Pr(p)(l - e-l & l‘)P 


which gives (54), by noting that 


Lemma A.2. If a > tr 2 /2 and b > 0, ffiere exists a constant C depending on a, b, a and T, such that 

.2 _2 rt 


sup Eexp f <C{ 1 + x 2 ), 

te[o,r] V 8 Jo r s) 


□ 


(55) 


where u = ^ — 1 > 0. 

<7^ — 
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v 2 a z f ds v 2 f ds 

Proof. For any t G [0,T], we set H t = %y/rf, so that Eexp(—-— / —) = Eexp( — / —j). The process 

8 Jo r s 2 J o H s 

{H t ,t G [0, T}) solves 

dH t = ~ 7 ^ tt - TjHtdt + dW t , H 0 = -\fx. 

\<j a 2J h t 2 a 

For any t G [0, T], we set B t = H t — H 0 — j * (|§ — i) -^4. Let (Z t . t G [0, J 1 ]) defined by 

/ ft b l 2 /*t \ 

= exp ( — / -H s dB s - — / iT 2 dsj . 


By the Girsanov Theorem, under the probability Q such that ^ 
Indeed (7T t ,i G [0,T]) solves 


dH t = 


cr 2 2 H 


= -g-, ( B tl t G [O,! 1 ]) is a Brownian motion. 

+ dB t: t <T, Hq = —yfx 


J 7 , 

2 a 1 \ dt 2 


and under Q, we note that (iT t ) is a Bessel process with index v = =% — 1. Moreover, by the integration by parts 
formula, f* 2H s dB s = H 2 — Hq — and 


b r r2 b 


2 pt 


ba 


Z t = exp -~H f — — H 2 ds H—+ t— < exp —=-x + T—x . 


10 

Now, denoting by E*® the expectation relative to ! 

- 2 rt ds 


, ba 




,2 pt 


expl T 


ds 

m 


. b ba 

< exp ( —x + T -2 


,2 rt 


exp ly 


ds 

m 


Let E^ denotes the expectation relative to P 2 ^-, the law on 


p+ td>+) 0 f ^ B esse i process with index v, starting 
at -y/x. The next step uses the following change of probability measure, for u > 0 (see Proposition 2.4 in [11]). 

,2 P 




cr(R s ,s<t ) 


(oR t \ u ( v 2 f l ds 

= lvij exp l' 


j(0) 

2 Jo RV ^ 


a(R s ,s<t ) 


where (R t , t > 0) denotes the canonical process on C(M , R + ). Then, we obtain that 


Eexp 


v 2 f ds 

~2~ ,/n 77! 


< exp [ —x + T^ ) E 

< exp ( -L -x + T^r ) E 


rH 

y/x 

( 0 ) 

\y/z 


exp 


v 2 f ds 

~2 Jo #!/J 


(oRtY' 
V2 yfr) . 


It remains to compute E 


" ( o"R \ 

f ^/— J ' 1 ^ — 0) be a two dimensional Brownian motion. Then 


E 


( 0 ) 


(<jRt V 


2 y/x 


E 


( W * 1 ) 2 + ( W . 2 + ^) 2 


' / \ 

and an easy computation shows that E^^- 

K v§) I ~ c ^ ^ + x 2 ). 


□ 
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B Proofs of Propositions 3.4 and 4.4 

Proof of Proposition 3.4. To simplify the presentation, we consider only the case when k(x) and h(x) are nil. For any 
e > 0 and x > 0, we define for all t £ [0, T\, the process — Xf), satisfying 

jr = i + T ^j^ds + [* r s j*' e dw s , 


with <pl = 6' (X s x + Oe J?’ £ ) dd and i/j e s = f* ^ 


add 


X*+e9J* 


Under (H3), the trajectories (Xf,0 < t < T) are 


strictly positive a.s. (see Remark 2.2). By Lemma A.l, f ip e s dW s is a martingale. Then Jf’ e is explicitly given by 

Jt’ e = ex P (^J o fads + f ( s dW s ( i’lfds'j . 

that fo 57 fxf dWs = 5 lo § (?) - So h^ds ; 


We remark 


J?’ e < CV —exp I - 


1 b_(X_s). 

2 X? 


and 


ds-x Os) ds + / N’s - 


2\/*0 




We upper-bound the moments E( Jf’ e ) a , a > 0. As b(x) > 6(0) — Kx. for any p. 






xexp| /oT : "^i 

and by the Holder Inequality for p > 1, we have 




E(Jf’ £ ) Q <C { E 


I 1 ) 2 ( P - 1 ) 


exp 


ctp 


2(P-1) 




/o 


X* 


f ° 2 *1 

)1 

\ 

Jo 4X? J 

) \ 

1 


Xf \ 2lf=TJ 


Then, for any 0 < a < 4, for any p > 1 such that op < 4, 

E( Jf’ e ) a < C \ E 


The same computation shows that for the same couple (p, a) and for any 0 < P < 


E 


_OT e ) c 

(*?) 


?lf +/3 I — 


X 2 l L J J 


(56) 


(57) 


which is bounded according to Lemma 3.1. Hence, by (56) for (a,p) = (2, 2), there exists a positive constant C such 
that 

E(X t x+e - X?) 2 < Ce 2 , V t £ [0, T], 

and Xf +e tends Xf in probability. We consider now the process ( Jf , t £ [0, T]) solution of (19). Applying the 

integration by parts formula in (20), we obtain that Jf = Jexp(f*(b'(X f) — + ^-^)ds), from which by 

using (H3), we have 


Jf < 



exp ( - [ (6(0) - ?)t^? ) exp(XT) < C 


4 '2X, 



(58) 
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Moreover, 


J x t - Jt’ e = f b'(XZ)(JZ - J^)ds + J* - J* a ' € )dW t 


2v^F 

+ jf ( 6 '(Xf) - 0*) J?*ds + jf (^= - 1>t)J5' e dW a . 

We study the convergence of E(Jf — J t x,e ) 2 as e tends to 0. We set £f := | Jf — Jf’ e |. By the Ito’s formula, 

E(ff) 2 =E f 2 b'(Xf)(£Z) 2 ds + E [ 2(b'(X a ) - - Jf' e )ds 

Jo Jo 


-E 


jf( 


V 2 ^ 


(Jf - Jf e ) + ( 


2y/Xf 


- 


Js. 


We upper-bound the third term in the right-hand side of the expression above: as Z _. > ip e a and ( 

2\/X * 2\/X 


c/y/xf. 


E (—£= - r s ) J*’ e < CE (—£= - </>«) 


k 2^xF 


(Jf’ e ) 2 N j 
' 2\/Xf Y ' S ' A /*F/ 


•V'l) < 


An easy computation shows that -i/Xf (— Z= — ipl) < i/e 


/X? 


Then, 


E ( 


2v^F 


-^K X,£ < C'-\/iE 




V (*?) ! 

2>5>4 




<Cv^, 


where we have applied (57) with (a,p, /3) = (|, |, ^ < 2— = |). By using the same arguments with (58), 


E 


V 2 ^F 


{ J x s - ■£•')( 


2y/X| 


- J s x ’ £ < E 




< 


(«(*© 


-E 








<CVc 


where we have applied (57) with (a,p,/3) = (|, \ = |). An easy computation shows that | 6 '(Xf) — <jf s \ < 

eJf ,e || 6 ,, || 00 . Coming back to the upper-bound of E(£f ) 2 , we have 


E(£f ) 2 <cj^ E(£f) 2 ds + Cv/ef + E Qf ^(£f) 2 Js 


(59) 


To conclude on the convergence, as e tends to 0, we use the stochastic time change technique introduced in [3] to 
analyze the strong rate of convergence. For any A > 0, we define the stopping time t\ as 


t\ = inf{s e [0 ,T\, 7 (s) > A} with 7 (t) = / 

Jo 

Then, by using the Lemma 3.1 with the Markov Inequality, 


f cr 2 (is 

4X? 


and inf 0 = T. 


P(t a <T) = P( 7 (T) > A) < exp(-^)E ^ex P < Cexp(~). 
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Choosing A = — log(e r ) for a given r > 0, we have that P(ta < T) < Ce 2 and 

E(£*) 2 <E(£*J 2 + Ce*. 

With (59), we can easily check that for any bounded stopping time r < T, 


and for t\, 


E(£f) 2 < exp {C(T - a)) |e QT ^(^) 2 ) + 

E(£f x ) 2 < C\E + C 0 yfe, 


for some positive constants Co and C\, depending on T. After the change of time u = 7 (s), we can apply the Gronwall 
Lemma 

E(£*J 2 < CiE + Co\fk < TC 0 v / ^exp(C'iA). 

With the choice r = (4Ci) -1 and A = — log(e r '), we get E(£f x ) 2 < TCq€X As T is arbitrary in the preceding 
reasoning, we conclude that E| Jf — Jf’ e \ tends to 0 with e for all t £ [0, T\. Consider now 

9{ X t ) 9{ X t ) / / \rx\ rx jx,e f J ( srx i TX,e\ j ix J(srx\ 

- 9 ( x t ) J t = / 9 \ x t + eaJ t ) da ~ J t9 K x t ) 

e Jo 

= ( Jt’ e Jt) I' g'( x t + eaJnda + J? f (g'(X? + eaJf^) - g'{Xf)) da 
Jo Jo 

:= A e + B e . 

EA e < hW? — Jf’ e |, which tends to zero with e. B e is a uniformly integrable sequence, g' is a continuous 
function. By the Lebesgue Theorem, as X : f ' :+f tends Xf in probability, IV tends to 0 with e. As a consequence, 

E(5(^J— )~9( x t ) ) tends tQ E [ g f X ?)Jf] when e tends to 0. □ 

Proof of Proposition 4.4. The proof is very similar to the proof of Proposition 3.4. Again, we consider only the case 
when h(x) and k(x) are nil. Let Jf ,e = -(Xf +e — Xf), given also by 


Jt’ = exp ( / fids + j i/ s d.W s - 


*j£W 


ds , 


with f e s = fg bfXf + 0eJf’ e )d9 and iff = f () l • F° r an A ^ function 9{ x ) with bounded derivative. 


we have 


g (xr e )-g(xn 


- g'(X?)J? = J x ^ f 1 g\Xf + eOJ^dd - 
J 0 

= (Jt’ e - Jt) t g\ x t + eOJndd + J* f 1 (g'(X? + edjn - g'(Xf)) dd 
Jo Jo 


:= A e + B e . 


E(Jf’ e ) a < exp(a||h , || 00 f)Eexp(a / ipldW s -- / {f e s ) 2 ds) and, by using Lemma 4. 2 (h), one easily concludes 

Jo 2 J o 

that E( Jf ,e )“ < C and consequently that Xf +f converges to X'f in ififl). Then, by applying the Lebesgue Theorem, 
E|B e | tends to 0. Moreover, E|A e | < \\g’\\ x '^E\Jf^^^Jf^. We can proceed as in the proof of the Proposition 3.4, 
to show that E| Jf’ e — J t x | 2 tends to 0, but now the moments E(Jf’ e ) a , a > 0 are bounded and theLemma4.1 ensures 
that the E\Xf\~ p ,p > 0 are all bounded. □ 
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C End of the proof of Proposition 3.2 

To compute (t, x), we need first to avoid the appearance of Jf (1) in the expression of (t, x). We transform the 
expression of j^(t, x) in (25), in order to obtain x) as a sum of terms of the form 


r T-t 


E J^exp ^^ P(Xf(l))dsj r(X%._ t (l))J%._ t (l) J 

+ £ * E | exp QT p(Xf(l ))du^ Jf(l)A(Xf(l)) J ds 


for some functions (3(x), T(x), A(*). In this first step, to simplify the writing, we write Xf instead of X® (1). Two 
terms are not of this form in (25): 


I = 2E |exp ^2 £ \'(Xf)ds j f(X£_ t )£ * b"(Xf) J*dsj 

II = 2E{ £ ‘ exp (2£ b'(Xf)d^j £(t + s,Xf)b"(Xf) QT b"(Xf)Jfd^j ds}. 
The integration by parts formula gives immediately that 


r T-t 


pT—t 


II = 2E < 


b"(xf)j: 


du 

dx 


(t + u, Xf) exp 2 b'(Xg)d6 b"(Xf)du ds 


By using again the Markov property and the time homogeneity of the process (Xf ), 

pT—t—s 


E 


exp [2 £ * b'(Xg)d6j f"(X$_ t ) j T s = E exp (2 £ * f"(X y T _ t _ s ) 


and, by using (24), 

rT—t 


I =2 


J 4 E £(Xf)Jf exp (2 £ V(X§)<W) + s, x£ ds 

~ 2 £ ‘ E { b"(Xf)Jf exp (2 £ b'(Xg)d£ 

£-(t + s + u, XI) exp (^2 £ b'(X y )dO^ b"(Xl) 


pT—t—s 


Conversely, 


l 


T—t—s 


E 


E 


du 

dx 




f T-t 

/ E 

J s 


du 

dx 


(t + s + u, X«) exp (2 £ b'(X y )ddj b"(X y ) 

(t + u, Xf) exp (2 J b'(Xf)de) b"(Xf) I T t 


du 


y=xg 

du 


and then 


I =2 


v=x* 


du > ds. 


I E £'(Xf)Jf exp (2 £ If (X$)<wj £±( t + s,x£ds 
- 2E |^ T ‘ b"(Xf)Jf i^£ ‘ £-(t + u, Xf) exp ^2 £ b' (Xg)d6^j b" (Xf)d^ ds j . 
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Finally, replacing I and II in (25), we get 


E | exp ( 2 ^ >Jtvx 


b'(X*)ds 


) f (3) W 






+ 1^ E {exp {2 jf b\X*)dv^j J* 


d^u du 

x(3-^(t + s, x *)b"(X:) + -(t + s,Xl 


ds. 


To eliminate J x , we introduce the probability Q 3 / 2 such that 


dP 


r Then 


=E 3 / 2 |exp (2 / T \'(x:)ds ) f^HxZ-JZ^JZ-t 


+ ^ E 3 / 2 {exp {2 b\X x u )di 


u 


X ( ^9^2 "I" S ’ )b"{X x ) + + S ’ 


ds. 


Again, we note that z[ 1 ’ 2 ' > Jf = exp ^ f* b'(Xf)duj 


and 


=E 3 / 2 | exp ' b'(Xf)ds^j / ( 3 ) (*r- t ) 

+ £ 1 E 3 / 2 {exp {3 J* b\X x )di 


Pit'll f)ni 

x (3 —(t + s,X*)b"(X x ) + —(t + s,Xl 


ds 


where we write X x instead of Jf x (l).. Finally, as C ^ (X x (l)) = £^(Jf x (|)), we obtain the following expression 

for 0M): 


d 3 u 

fe 3 


r T-t 


(t, x) =E < exp 3 


tf(X?(|))d S ]/ (3) (Xf_ t (|)) 


J q T *E{exp (zj\\X* u (\)d^ ^^(t + s,Xf(^))b"(Xf^)) 


+^(t + s,Xf^))b^(Xf^)))\ds. 
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J®(|) exists and is given by (20). By the Proposition 3.4, x) is continuously differentiable and 


=E I exp | 3 


r T-t 


b'(Xg( — ))ds 


X 
c T-t 


+ [ ‘sjexp (3^{t + a,X^W{X^)) 

+ g(i + s, jT 36 W (X;(|)) JJ(|)d U } ds 

® |exp (z j\'{X* u {\))d^ r s {\) (s0(t+ S ,^(|))&"(^(|)) 

+ g( i + S ,X-(|))&( 4 )(X-(|)))|d Sj 


from which we can conclude on (16). 
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